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1. Introduction 

In hadronic collisions, jet cross sections are computed as a convolution of hard par- 
tonic cross sections with parton-distribution functions, followed by the fragmentation 
of the final-state partons into hadrons. The theoretical prediction can be improved 
by including higher order corrections which have the effect of reducing the unphysical 
renormalisation- and factorisation-scale dependences and by improving the matching 
of the parton-level theoretical jet algorithm with the hadron-level experimental jet 
algorithm. At present jet production is described at next-to-leading order (NLO) 
and several numerical programs are available [1, 2] which have been extensively used 
to compare with data from the TEVATRON and CERN SppS. 

Improving the theoretical prediction to next-to-next-to-leading order (NNLO) 
requires several ingredients. First, the parton-density functions are needed to NNLO 
accuracy which in turn requires knowledge of the three-loop splitting functions. At 
present, the even moments of the splitting functions are known for the flavour sin- 
glet and non-singlet structure functions F 2 and F L up to N = 12 while the odd 
moments up to N = 13 are known for F% [3, 4]. The numerically small Np non- 
singlet contribution is also known [5]. These moments are sufficient to parameterise 
the splitting functions in x-space [6, 7] and NNLO global analyses [8] are starting to 
appear. Second, the hard scattering matrix elements can be computed at NNLO. In 
the high energy limit, the quarks can be assumed to be massless and at this order, 
there are contributions from 2-^4 tree-level diagrams [9, 10, 11, 12], from 2^3 
one-loop- level diagrams [13, 14, 15] and two- loop 2 — > 2 diagrams. 

The evaluation of the two-loop diagrams has been a challenge for the past few 
years due to the presence of two-loop planar and crossed boxes. In the massless 
parton limit and in dimensional regularisation, analytic expressions for these basic 
scalar integrals have now been provided by Smirnov [16] and Tausk [17] as series 
in e = (4 — D)/2, where D is the space-time dimension, together with constructive 
procedures for reducing tensor integral to a basis set of known scalar (master) in- 
tegrals [18, 19]. This makes the calculation of the two-loop amplitudes for 2 — > 2 
QCD scattering processes possible. Bern, Dixon and Kosower [20] were the first to 
address such scattering processes and provided analytic expressions for the maximal- 
helicity-violating two-loop amplitude for gg — > gg. Subsequently, Bern, Dixon and 
Ghinculov [21] completed the two-loop calculation of physical 2—^2 scattering am- 
plitudes for the QED processes e + e~ — > fi + fi~ and e + e~ — > e~e + . More recently, 
we have provided expressions relevant for unlike- and like-quark scattering in the 
massless limit in Refs. [22, 23] respectively. The corresponding matrix elements for 
the self- interference of one-loop 2^2 quark scattering processes are given in [24]. 

In this paper, we address the O (af) one- and two-loop corrections to the QCD 
process 

q + q^g + g, (1.1) 
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together with the time-reversed and crossed processes 



q + g -> q + g, 

g + q -> g + q, 
g + g -> q + q. 



(1.2) 
(1.3) 
(1.4) 



As is in Refs. [22, 23, 24], we use the MS renormalisation scheme to remove the ul- 
traviolet singularities and conventional dimensional regularisation, where all external 
particles are treated in D dimensions. We provide expressions for both the inter- 
ference of tree-level and two-loop graphs as well as the self-interference of one-loop 
amplitudes. In each case, we find that the infrared pole structure agrees with that 
obtained using Catani's general factorisation formulae [25]. The finite remainders 
are the main new results presented in this paper and we give explicit analytic ex- 
pressions valid for each of the processes of Eqs. (1.1)— (1.4) in terms of logarithms and 
polylogarithms that are real in the physical domain. For simplicity, we decompose 
our results according to the powers of the number of colours N and the number of 
light-quark flavours Np ■ 

Our paper is organised as follows. We first establish our notation in Sec. 2. 
Analytic expressions for the interference of the two-loop and tree-level amplitudes 
are given in Sec. 3, while formulae describing the self-interference of one-loop graphs 
are given in Sec. 4. In Sec. 3.1 we adopt the notation used in Ref. [25], to isolate 
the infrared singularity structure of the two-loop amplitudes in the MS scheme and 
we demonstrate that the anticipated singularity structure agrees with our explicit 
calculation. The finite O (e°) remainder of the two-loop graphs is one of the main 
results of our paper and expressions appropriate for the qq — > gg, qg — > qg (gq — > gq) 
and gg — > qq scattering processes are given in Sec. 3.2 in terms of logarithms and 
polylogarithms that have no imaginary parts. Expressions for the self-interference 
of one-loop graphs are given in Sec. 4 in terms of the one-loop bubble integral in 
D = 4 — 2e and the one- loop box integral in D = 6 — 2e. Analytic formulae connecting 
these integrals in the various kinematic regions are given in Appendix A. As for the 
two-loop contributions, the singularity structure agrees with that obtained using 
Catani's formalism. Finally we conclude with a brief summary of the results in 
Sec. 5. 

2. Notation 

For calculational purposes, the process we consider is 



where the particles are all incoming and carry light-like momenta, satisfying 




(2.1) 



Pi 



+p%+p%+p i Z = o, p-=o. 
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The associated Mandelstam variables are given by 

S=(P1+P2) 2 , t=(p 2 +p 3 ) 2 , U=(p 1 +p 3 ) 2 , 



s + t + u = 0. (2.2) 



We work in conventional dimensional regularisation treating all external quark and 
gluon states in D dimensions and renormalise the ultraviolet divergences in the MS 
scheme. The bare coupling a is related to the running coupling a s = a s (fi 2 ) at 
renormalisation scale /i, by 



(2.3) 



where 



(4tt)V 



£7 



7 = 0.5772 . . . = Euler constant 



(2.4) 



is the typical phase-space volume factor in D = 4 — 2e dimensions. The first two 
coefficients of the QCD beta function, (3q and for N F (massless) quark flavours, 
are 

\1C\ - 10C A T R N F - QC F T R N F 



IIC A -AT R N F 
6 
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where N is the number of colours and 

N 2 -l 



2N 



C A = N, 



T, 



R 



1 

2' 



(2.5) 



(2.6) 



The renormalised four point amplitude in the MS scheme is thus 



\M) = 47ia s 



\M®) 



2vr 



|A^ (1) ) 



2tt 



\M (2) } + 0(af) 



(2-7) 



where the \M^) represents the colour-space vector describing the i-loop amplitude. 
The dependence on both renormalisation scale \i and renormalisation scheme is im- 
plicit. 

We denote the squared amplitude summed over spins and colours by 



(M\M) = J2\M(q + q^g + g)\ 2 = C(s, t, u). 



(2.J 



which is symmetric under the exchange of t and u. 

The squared matrix elements for the crossed processes are obtained by exchang- 
ing the Mandelstam variables and introducing a minus sign for each quark change 
between initial and final states 



J2\M(q + g 



q + q)\ 2 
i + g)\ 2 
g + q)\ 2 



C(s,t,u), 
-C(u,t,s), 
—C(u, t, s). 



(2.9) 
(2.10) 
(2.11) 
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The function C can be expanded perturbatively to yield 



C(s,t,u) = 16ttV 



C 4 (s, t,u)+( C 6 (s, t, u) + (^) 2 C s (s, t,u) + (a?) 



2tt7 ^ ' ' ' V2tt 



(2.12) 



where 



C 4 (s,t,w) = (M (0) \M {0) ) 



N 2 -l, , / iV 2 - 1 N' 
= 2 — — — 1 - e 



-2^ (t 2 + M 2 -e S 2 ), (2.13) 



N ' \ ut 

C 6 (s,t,u) = ((M (0) |M (1) > + (M (1) \M {0) )) , (2.14) 
C 8 (s,t,u) = ((M (1) \M (1) ) + (M (0) \M^) + (M (2) |M (0) >) . (2.15) 

Expressions for C 6 are given in Ref. [26] using dimensional regularisation to isolate 
the infrared and ultraviolet singularities. 

In the following sections, we present expressions for the infrared singular and 
finite contributions to C 8 and the crossed processes. For convenience, we divide 
C 8 (s,t,u) into two pieces 

- the pure two-loop contributions 

C 8 ^°\s,t,u) = (MW\M<V) + (M<V\M<®), (2.16) 
described in Sec. 3 and 

- the self-interference of the one-loop amplitude 

C s( - lxV >(s,t,u) = (M^\M m ), (2.17) 

described in Sec. 4. 

As in Refs. [22, 23], we use QGRAF [27] to produce the two-loop Feynman dia- 
grams to construct \M^). We then project by (M^\ and perform the summation 
over colours and spins. It should be noted that when summing over the gluon po- 
larisations, we ensure that the polarisations states are transversal (i.e. physical) by 
using an axial gauge 

V ete? = -g^ + n *rt + n irf (2.18) 

" Tls • Pi 

spins ' r t 

where pi is the momentum of gluon % and rij is an arbitrary light-like 4-vector. For 
simplicity, we choose = p% and = p%. Finally, the trace over the Dirac matrices 
is carried out in D dimensions using conventional dimensional regularisation. It is 
then straightforward to identify the scalar and tensor integrals present and replace 
them with combinations of the basis set of master integrals using the tensor reduction 
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of two-loop integrals described in [18, 19, 28], based on integration-by-parts [29] and 
Lorentz invariance [30] identities. The final result is a combination of master integrals 
in D = 4 — 2e for which the expansions around e = are given in [16, 17, 18, 19, 28, 
31, 32, 33, 34]. 

3. Two-loop contribution 

We further decompose the two-loop contributions as a sum of two terms 

C 8(2x0) (s,t,w) = Poles(s,t,u) + Finite(s,t,u). (3.1) 

Voles contains infrared singularities that will be analytically canceled by the in- 
frared singularities occurring in radiative processes of the same order (ultraviolet 
divergences are removed by renormalisation). Finite is the remainder which is finite 
as e — > 0. 

3.1 Infrared Pole Structure 

Following the procedure outlined in Ref. [25], we can write the infrared pole structure, 
renormalised in the MS scheme as 



Voles(s, t, u) — Re j 

+ — ( — A(e, s, t, u)B(e, s, t, u) — B(e, s, t, u)C(e, s, t, u) + 21ZB(2e, s, t, u) 



4— B(e, s, t, u)) + j^B(e, s, t, u)D(e, s, u, t) 



V 



AN 



%.{s,t,u) 



+ 



lA(e, s, t, uf + B(e, s, t, u)D(e, s, t, u) + 4yA(e, s, t, u) - 2TZA(2e, s, t, u) 



+ — ( — A(e, s, t, u)B(e, s, t, u) — B(e, s, t, u)C(e, s, t, u) — 4— S(e, s, t, 



u 



+ 
+ 
+ 



\ V 2 

+21lB(2e, s, t, u)j - J^ B (^ s > t, u)D(e, s, u, t) 

V 2 V l 

— A(e, s, t, u) + jB(e, s, t, u) d{s, t, u) 

A(e, s, t, u) + jB(e, s, t, «) £ 2 (s, t, u) 

V 2 V V 

— D(e, s, t, u) + — C(e, s, t, u) - ^^(e, s, u, f) 



T 2 (s,t,u) 



C 3 (s,t,u) + H 2 



+ (u<->M, 



(3.2) 
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where 



A(e, s, t,u) = — 



+ N 



B(e, s, t, u) 
C(e, s, t, u) 
D(e, s, t, u) 



1 3 
+ — + 



K 2Ne 4e j 



e 2 4e 2Ne 



u 



.01 
t 



+ c F - 



A 4 



2\ e 



+ 



i/. 



2\ e ' 



^7 



r(i-e) 



.01 

s 



=«7 



2e 



iVe 



r(i-c) 

2 \ e P £7 



s ; r(i -e) 



and 



e 2 4e 2iVe 



V = AT 2 - 1, 



.01 
t 



.01 

u 



3 «7 



r(i- e ) 



r(l-2e) fa 



r(l-e) V e 



X = 



67 7T 2 



10. 



18 6 / ° A ~ ~9 TrNf - 



The tree-type structures, 7[ and 7^, are given by 

7i(s,*,«) = -T 2 (s,t,u), 



u 



T 2 (s,t,u) = 8(1 -e 



(t 2 + u 2 - es 2 ) 



(3.3) 
(3.4) 
(3.5) 
(3.6) 

(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 



while the interference of tree with one-loop structures are represented by Ci, C 2 and 
N 2 + 1 



d(s,t,u) 



v L /2(S, t, U) - M > + fs(s, t, U) 

-6A) 77 ~^ Ti(s,t,«) Bub(s) 



3-2e 



-2e(l-2e) 



AT 



: Bub(«) + 



iV / 1_ _ 2^o 3 
2 I e 2 Ne 2e 



Bub(s) 



1 



2iV Ve 2 2e 
C 2 (s,t,u) = -Ci(s,u,t) 



+ — Bub(s) 



Ti(s,t,w) 



£ 3 (s, i, u) = 16(1 - 2e)- [t 2 + w 2 + (u* - 2t 2 - 2m 2 ) e + (t 2 + m 2 + 3uf) e 5 
+2 f- -2) [Bub(f) -Bub(s)]7!(s,t,u) + /i(s,t,u) 



(3.12) 
(3.13) 
Box 6 (t,w) 



+2 f- - 2) [Bub(u) - Bub(s)] T 2 (s,t,u) + ( 3 - 14 ) 
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where the infrared- finite functions /i, fa and / 3 are 

fi(s,t,u) = 8(1 - 2e)^[u (2u 2 + 5t 2 + 3tu) + e (3t 3 - 4m 3 - 3tu 2 ) 

-e 2 s (At 2 + 2m 2 + htu) + s 2 te 3 ]Box 6 (s, t) (3.15) 
f 2 (s, t, u) = 8(1 - 2e)^ [(t - m) (t 2 + 2m 2 + tu) + e (-2t 3 + Mt 2 + Atu 2 + 5m 3 ) 

(3.16) 



-s 3 e 2 



us' 
Box 6 (s,m) 



f 3 (s,t,u) = -4(1 -e) 



sm liV 



2s - t - e(2s - 3m) - 3se 5 



+4iV(w - ue + se 2 )} [Bub(u) - Bub(s)] 



At 



[-N [18m 2 + 15t 2 - 3t(s - t) 



s 2 w(l-e)(3-2e) 
-e (78m 2 - 36t(s -t) + st) + e 2 (80m 2 + 10s(s - t) - 69st) 
1 



N 



-24m 2 + 3tu - 21t 2 + e (85m 2 - 43t(s - t) + 3st) 



20s 2 - 40tw 



-e 2 (112m 2 + 6s(s - t) - 109st)] + (3 

-2e (38s 2 - 3ms - 62tw) + 4e 2 (27s 2 - 26fu)] } eBub(s) + C(e 3 ). 

(3.17) 

These expressions are valid in all kinematic regions. However, to evaluate the pole 
structure in a particular region, the one-loop bubble graph Bub and the one-loop 
box integral in D = 6 — 2e dimensions, Box 6 , must be expanded as a series in e. This 
analytic expansion is given in Appendix A. 

The function 7i 2 , that appears in Eq. (3.2), exhibits only a single pole in e and 
is given by 



7i 2 ^(M^\H {2 \e)\M^) 

where the constant is 

^ (2) = (C3 + ^ + ^vr 2 )cl+(l3C 3 + 



^7 



4er(l - e 



245 23 



(3.18) 



7T 



58 



108 24 

7T 2 29 N 



--7T 2 C A C F + 



10 
27 



Nl 



+ " 7^ )C A N F + — + — iV F C F + - 7 + 7T 2 - 12 C 3 )C F , 



36 27 



12 54 



(3.19) 



and ( n is the Riemann Zeta function with ( 2 = 7r 2 /6 and ( 3 = 1.202056 . . . We note 
that is renormalisation-scheme dependent and Eq. (3.19) is valid in the MS 
scheme. We also note that Eq. (3.19) differs from the corresponding expressions 
found in the singularity structure of two-loop quark-quark scattering in all but the 
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Cp coefficient. This is due to the presence of infrared emissions from gluons which 
modify the terms involving either Ca or Np. 

It can be easily noted that the leading infrared singularity in Eq. (3.2) is O (1/e 4 ). 
It is a very stringent check on the reliability of our calculation that the pole structure 
obtained by computing the Feynman diagrams directly and introducing series expan- 
sions in e for the scalar master integrals agrees with Eq. (3.2) through to O (1/e). 
We therefore construct the finite remainder by subtracting Eq. (3.2) from the full 
result. 

3.2 Finite contributions 

The finite two-loop contribution to C 8 (s,t,u) is defined as 



In hadronic collisions, all parton scattering processes (Eqs. (1.1)— (1.4)) contribute 
simultaneously. We therefore need to evaluate Tinite{s ) t, u) for the qq — > gg and 
99 QQ process (which we denote as the s-channel since, although the tree-level 
process contains graphs in all three channels, the squared tree matrix elements are 
proportional to 1/s 2 ) and Finite{u,t,s) for the QCD Compton processes qg — > qg 
and gq — > gq (which we label as the w-channel) . 

Of course, the analytic expressions for the various processes are related by cross- 
ing symmetry. However, the master crossed boxes have cuts in all three channels 
yielding complex parts in all physical regions. The analytic continuation is therefore 
rather involved and prone to error. We therefore choose to give expressions describ- 
ing C 8 (s,t,u) and C 8 (u,t,s) which are directly valid in the physical region s > 
and u, t < 0, and are given in terms of logarithms and polylogarithms that have no 
imaginary parts. 

As usual, the polylogarithms \A n (w) are defined by 



Using the standard polylogarithm identities [35], we retain the polylogarithms with 
arguments x, 1 — x and (x — l)/x, where 




(3.20) 




(3.21) 



_ ! 

s 




u 



X 



- 1 



(3.22) 



y 




X 



For convenience, we also introduce the following logarithms 




(3.23) 
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where /x is the renormalisation scale. 

For each channel, we choose to present our results by grouping terms according 
to the power of the number of colours N and the number of light quarks Np, so that 
in the generic c-channel we write 



Finite c (s,t,u) = (iV 2 — 

x (n 3 A c + NB C + lc c + -^D c + N F N 2 E C + N F F C + ^G c + N F NH C + . 

(3.24) 

3.2.1 Finite(s,t,u): the s-channel process 

We first give expressions valid for the annihilation processes, qq — > gg and gg — > qq. 
We find that 

{( 36077 109 4 241 5 2 2327 „ . , 11 , „ n ^ „ 

^ = {(-^"144' + l2<*-2* 2 ^S + 20Lu( X ) + -^S-9( 3 S 

+ X -f S 2 - 21 Li 4 (y) - 20 Li 4 (z) - ^ Li 3 (y) + ± X F 3 + g F 2 X + H y X vr 2 

+4 7T 2 Li 2 (x) + 16 F ( 3 + ^ F 2 Li 2 (y) + ^ F Li 2 (y) + ^ F 5 - 5 X 2 Y 2 

+ ^ Y 2 -5X 2 --Y 3 + -Y 4 --X 4 + -YX 3 --7T 2 Y 2 --n 2 X 2 
72 36 8 6 3 3 3 

_ii! y v-Ii>-s-20YLi 3W + ^rU 

36 3 216 y w 

/ll 4 9, 169 2 457 11 „ 33 VT .,, n „ > 33 aa2 „ 

-f l »m - w y + 1 y * 2 + 8 y Ll » w " s - Y3 + if y2 + f y s + - y2 y2 ! 

/29 4 205 A o 2 1535 30377 nrtT ., , 11 2 „ „ n , 

949 ^^79^ 
• S 2 + 18 Li 4 (2/) + 16 U 4 (z) + Li 3 (y) + ^ Li 3 (x) - - X F 3 - F 2 X 



9 " v ' 12 w/ 4 v ' 3 24 

-| F X vr 2 + ^ X - 4 vr 2 Li 2 (x) - F 2 Li 2 (y) - A y Li 2 (y) - | X 5 

+ | 5 X 2 - 8 Li 3 (y) X - F 5 + 2 X 2 F 2 - ^ X Li 2 (x) - " X 2 F 
+ 4XL l3 (,) + 4X C 3-X 2 L l2 (,) + fx. 2 -^F 2 -!gx 2 + ^y3 

-iF 4 + lx 4 + -^X 3 -^FX 3 + ^ 2 F 2 + 4 7 r 2 X 2 + ^F 7 r 2 + |F 2 5 



+8FLi 3 (x)--^F)^ + ^^^ (3.25) 
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{ [Wo " 4 - ll Cs " '-W " 2 " W 5 + 2 F Lhiy) + 8 Li4(x) ^ 20 F X Li2 (y) 

y X S + ^ vr 2 S + 10 C 3 S - ±^ S 2 - 79 Lu(y) - 74 Li 4 (z) - ^ Li 3 (y) 
25 T ... 164771 5 v - 1 . 19 ^ v 2 22 ^ 19 2t . , . 

+ y Ll3( ^ + w + 3 xr + i r x+ y rX7r "27 x -y 7r Ll2(a;) 

+74 Y C 3 + ^ Y 2 Li 2 (y) + ^ F Li 2 (y) + ^ X F - ™ X S + 24 Li 3 (y) X 

- f V S - y X 2 F 2 - | X Li 2 (x) + y X 2 y + 44 X Li 3 (x) - 36 X ( 3 

-llX 2 L l2 (x)-^X. 2 -^F 2 + ^X 2 + ^F3_3 r4 _25 x4 _^ x 3 
v ; 36 72 36 36 8 12 2 

+9 V X 3 — — 7r 2 Y 2 — — 7r 2 X 2 — — y 7r 2 + — y 2 S — 70Y U 3 (x) -^y)- 
3 3 36 3 216 J u 

( 229 , 118 , 799 2 11 „ 139 in ^ T . /N 65 T . . , 

+ vr 4 + Cs vr 2 5+ 18 y Li 3 (y) + — UAy) 

\ 240 3 S 72 6 36 KyJ 2 Kyj 

9QQ 9 ^ 14 7 7S 

" Li 3 (i) + - X Y 3 - i Y X ir 2 + ^ X + 19 Y ( 3 + '- Y 1 Li 2 (y) +^-XY 



AyS-—X 2 Y 2 + ^-X Li 2 (x) + 5? A - Y-—Y'-- Y< + ™ A" 
6 8 6 V ; 3 36 48 18 

47 , ~ 55 _ 11 



"12 ^ 6 ^ 77 6 " 9 y Li3(;r) J 
+ ( " + ^ Cs + f vr 2 + ^ 5 + 20 y Li 3 (y) + I Li 4 (x) + 8 Y X Li 2 („) 

^ - 2 Cs ^ - y Li 4 (y) - 22 U 4 (z) + ^ Li 3 (y) - ^ Li 3 (x) -ylY 3 

+ f Y 2 X - | y X vr 2 + ^ X + 7 vr 2 Li 2 (*) - y Ca - f V 2 Li 2 (y) - ? j- Y Li 2 (y) 

+^XY+^XS + ^SX 2 + 7U 3 (y)X-^YS + ^X 2 Y 2 + 7 -±XU 2 (x) 

-24 X 2 y + 16 X Li 3 (x) + X C 3 - y X 2 Li 2 (x) + ^ X vr 2 - ^ y 2 - ^ X 2 
101 - 31 , 4 33 341 o ota „ 11 - „ 11 2 - 103,, 2 

y 3 + — y 4 — x 4 + — x 3 + 3yx 3 + — ?r 2 y 2 + — ?r 2 x 2 y?r 2 



18 48 16 18 4 4 18 



11, , 2 C — ' * 19139 601 ^ t<2 
t 4 



-— y 2 5- 7yLi 3 (a;) - + — Y I ^ 

6 V ; 324 36 / s 2 



- X2 ^) + \ u ^ t ) ( 3 - 26 ) 



41393 301 4 173 , 133 2 5297 „ ,„ T ., S ™ t -,n 



2592 720 36 36 216 
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13 15 
-4 Y X Li 2 (y) - — 7T 2 5 + 5 Cs S + Li 4 (y) + 8 Li 4 (z) - — Li 3 (y) - 3 Li 3 (x) 

-f* " ^ ^ X - f Y X n 2 + | X - jjvr 2 Li 2 (x) - 12 F (s - \y 2 Li 2 („) 

1 5 22 T 

+ — Y U 2 (y) + 2XY - — SX 2 + 10 Li 3 (y) X - 11 Y S + - X 2 Y 2 + 3 X Li 2 (x) 

10 , / \ , , s 49 9 425 , 179 , 

- - X 2 y + 10 X Li 3 (x) - 6 X C 3 - 4 X 2 Li 2 (x) + — X n 2 - — Y 2 + — X 2 

_H!y3 + 3 y4 + I X 4_ 65 X 3_8 yX 3 + 5 7r 2y2 + 3 2 X2 _5_ y 2 

36 8 3 18 3 6 2 18 

-H^ + 14 y Li3(l) + |,)i 

+ " 4 + " ¥ " 2 + T X Li2{x) + 26Y Lkiy) ~ T y2 Li2(y) " 25 F Cs 

119 13 T ... 3133^ 4 _~ 7 2 - 55^ 2 55 v ^ 

X 2 Y - — Li 3 (x) + —— X - - XY A + -it 2 Y 2 + — Y n 2 - —XY 



24 4 ° v ' 144 3 4 18 12 

+ l^ 4 -!^ + !H 



f 13 „4 , , X „2 3 C , 11 T . , , , 11 



- (— n i + 15( 3 + -7r 2 --S+-U 4 (x)+n 2 S-12( 3 S--U 4 (y)-6U 4 (z) 

+11 Li 3 (y) - - Li 3 (x) - ^- ^y 2 X + -yX7r 2 --X + l7r 2 Li 2 (x) 
4 vy; 4 v ; 16 8 3 16 3 v ; 

+3 y C 3 - \ Y 2 Li 2 (y) - ^ y Li 2 (y) - y X y + 3 Li 3 (y) X - 1 X 2 Y 2 

+ Li 2 (x) + f X 2 y - 3 X C 3 + ^ X 2 Li 2 (x) - I X vr 2 + g y 2 - 11 X 2 

-3yL l3 (,) + §y)J + X 2 -^U{^t} (3.27) 



{(-^ 7r4 + ^3 + ^ 2 -^ + 2yLi 3 (y) + \n 2 S-6C 3 S + llU 4 (y) 

3 1 87 3 7 

+14 Li4(z) - - Li 3 (y) + 8 Li 3 (x) - — + - Y 2 X - - Y X n 2 + X + n 2 Li 2 (x) 

- 18 y C 3 - ^ V- 2 Li 2 (y ) + ^ y Li 2 (y ) + 2 Li 3 (y ) X + 4 X 2 Y 2 - 8 X Li 2 (x) 

-11 X 2 y - 6 X Lia(x) + 4 X ( 3 - X 2 Li 2 (x) + ^ X vr 2 + | y 2 - \ X 2 - \ Y 3 
2 o o 4 4 

_ly4_1^4_ Z yx3 1 2y2 1 2x2 2 ( ^ 39 \ 

8 12 3 2 2 3V ; 8 / 



11 



+ (i - i£ 71-4 - 2 Cs + § 7r2 - 19 x Li2(;r) - 4 r Li3( ^ - y2 Li2(?/) - 12 F C3 

-^X*Y + 19U 3 (x) + ™X-^Y 2 + ^Yn*- 7 -XY + 10U 4 (y) 

-- Y X TV 2 + 16 Y U 3 (x) - — X 3 + 4 Y 2 - — Y A + 2 X 2 Y 2 ] 
3 V ; 6 12 J 

+ ( -3X + 7 -X 2 + 3Y + ^tt 2 + \y 2 -3XY] 



t 2 



s 2 



t 4 



-X 2 — + U~t (3.28) 



J (14 4 17 77 2 185 1 2 44 2 

< — n OH 7r H o — 8 Li 4 (x) 7r S 6+8 LiAy) 

\\45 6 S 54 54 v ; 36 9 KyJ 

+8 U 4 (z) + ^ Li 3 (y) - | F 2 X - 1 y X tt 2 - ^ tt 2 Li 2 (x) - 8 Y ( 3 - ^Y Li 2 (y) 

-y y 5 + 2 x 2 y 2 - ^ y 2 + 2 x 2 + ^ y 3 + i x 4 - ^ y x 3 + ^ 7T 2 x 2 + j| y 7T 2 

2 , r2 „ , , 10 _ 1307\ i 



+ -y< S + SYU,W + T Y + — 



u 



+ [-% - ~^^-^ 2 + \S -AXU 2 {x) + 2Y C,-2X 2 Y + AU,{x) 

+^X+ 2 -Yn 2 -2YU,{x) -2Y 2 - ^-Y S- \x 2 Y 2 ^ 

+ (-— 7r 4 + -C 3 --vr 2 -— S + 4U 4 (x) + -n 2 S+-S 2 -4U 4 (y) 
\ 180 3 S 27 27 v ; 18 9 KyJ 

4 2 2 19 2 

-4 U 4 (z) + | Lia(y) - 4 Li 3 (x) - | y 2 X + | Y X tt 2 - ±| X + | tt 2 Li 2 (x) 

+2yC3-^Li 2 (y)-^ + ^X5-^X 2 + 2Li 3 (y)X+ 1 ^y5 

-ix 2 y 2 + 4XLi 2 (x) + 2X 2 y-2XC 3 -— Xtt 2 + — y 2 + Ix 2 --y 3 

2 9 9 9 9 

_Ix 4 --x 3 + -yx 3 - -7r 2 x 2 --y7r 2 --y 2 5-2yLi 3 (x)--y^ 41 

6 9 3 3 9 3 v; 6/s 2 J 

+ Jw^ti (3.29) 



-16 Li 4 (y) - 16 Li 4 (z) + 1 Li 3 (y) + 4 Li 3 (x) + ^ Y 2 X + y X tt 2 - ^ X 



12 



+ l^U 2 ( X ) + l Q Yt 3 -lYU 2 (y)-lxY-^ + fxS + 2 {YS 

-4X 2 Y 2 -4XU 2 (x)--X 2 Y--Xn 2 + -Y 2 --X 2 --Y 3 --X i 
v ; 3 18 18 18 9 3 

+ 2 -X 3 + 8 -YX 3 -K 2 X 2 + ^Yn 2 -^Y 2 S-16YU 3 (x) + ^ Y ) *- 
33 3 18 3 9 J u 

+ — tt 4 + — C 3 + -7T 2 + -S +-XU 2 (x) +8Y( 3 + ^X 2 Y--U 3 (x) 

V 90 3 S 12 3 3 K J ^6 3 K J 

-^X + ^Yk 2 -\xY + \y 2 S-SYUz{x)-^X* + ^Y 2 - 1 4; 
9 18 3 3 w 36 18 36 

+^YS-X 2 Y 2 ^J 

J _^ S --XS+-Xn 2 --SX 2 + -X 2 Y + -YS--n 2 S 

V 27 3 2 3 2 36 

l v2n 7^ 4 t • / \ 211 - 4 , , 4085 4 T . /N 1 
"^3 9 3 L " W 3 L " W + W + 3 Ll3 W " 5 <' 

9 3 vy; 2 36 6 36 18 18 

-?^)?-^}+H} < 33 °) 

/ ^ + 143 W 2 _ 227 ^ 2 

\\15 18 S 36 54 v ; 12 Kyj v ; 

-8U 3 (x) + ^ + lY 2 X-8YXn 2 + ^X-8n 2 Li 2 (x)-48Y( 3 

+ ixr + -SX 2 + 2Y S + 12 X 2 Y 2 + 8XLi 2 (x) + ^X 2 F + — Xvr 2 
2 3 3 9 

+ - Y 2 + — X 2 + — Y 3 + 2 X 4 + - X 3 - 8 Y X 3 + 4 tt 2 X 2 + - Y ix 2 + - Y 2 S 
9 9 18 9 9 3 

+A8YU 3 (x)-^-Y ) - 
3 J u 

/11 107 71 71 

+ f - tt 4 - 52 C 3 - ^ tt 2 - 52 X Li 2 (x) - 24 Y C 3 - y X 2 Y + 52 Li 3 (*) + y X 

+ ^Y it 2 - \xY + \y 2 S + 24YU 3 (x) + \x 3 - ^-Y 2 + \y S + 3X 2 Y 2 \ 
y o o y y o / 

+ ^ - y X tt 2 - 8 X - 4 X Li 2 (x) - 18 X 2 + 4 Li 3 (x) - 4 Li 3 (y) + 8 Y + 4 F Li 2 (y) 

+ T y ' 2 + 9y2 )S + 9A ' 2 ^} + {"« i } < 3 - 31 ' 

I I 27 9 27 27 9 6 J u 



13 



.2 ,40 G 8 o2 , 20 „ 8 „ 1 A ( 2 



+ 27^ + 27 5 -9 5 + 27 y -9 yS -3 y b + (3 ' 32) 



\ V 9 27 6 27 9 9 18 54 

-5 y ' s +2H^H A -- y )S}+H} (333) 

3.2.2 J-inite(u,t, s): the w-channel process 

Similarly, for the Compton scattering processes, qg — > qg and — > gg, we find that 
the coefficients for the finite part in Eq. (3.24) are given by 

/ 2417^ 11 4 19 2 3 A 11 33 T .,. 457 OT . /w 

"f ^ y - f %2 Y + 8 Cs Y ~ 1 " 2 y2 ~ 2 A ' 2 y - T x hh(x) ~ 8 x 6 

_|,, 2 + i, y .3 + 8U3fe) , + ff/ , + -, 2 _iZ, 3 _lll y , 

8 2 4 141 2417 \ 

+2F 2 X + 2X 2 r 2 --7r 2 XF--F 4 --F 3 + — Y 2 + —Y) 

3 3 3 8 72 / 

+ (m 36077 121 2273 181 215 103 

1 9 864 ^ 72 432 720 144 6 S3 



+ If ^ + ^ Li4(z) + 1 UM + 1 Li3(a;) " 10 Lis {x) Y ~ l ~Y u Y ~ 1 " 2 r 

-6C3y + ^7r 2 y 2 + ^ 2 y-gxLi 2 (x) + 8XC3-^ 2 Li 2 (x)-|x7r 2 
-y C/ X 2 - H X F 3 + 1 vr 2 X 2 - \ vr 2 Li 2 (*) - y £/ F 2 - 10 Li 3 (y) X 

+ 1 ^ux + ^x* + 1 -^x 2 - 7 1x*- 1 -x*y- 1 -^yx- 7 -y 2 x- 1 1x 2 y 2 

18 12 144 72 2 72 4 8 

+ n" 2 * y + t u x Y ~\ ^ » x Y + i§ y 4 + 1 Y 3 + f y 2 



2273 > 



216 



t 2 + s 2 
si 



/ll rr 2 2273 v 35 4 865 2 35 , 41 . « 35 T . . , 

+ — [/ 7T 2 X 7T 4 7T 2 + — C 3 Li 4 (z Li 3 (y) + 20 Lii(x) 

V 6 432 36 144 12 s3 2 v ; 6 vy; V ; 

4135 T 13 

+ y Li 4 (y) - — Li 3 (x) - 10 Li 3 (x) Y - - vr 2 Y - 20 Li 3 (y) Y + 10 ( 3 Y - — vr 2 Y 2 

T o 35 , 1 o 55 o 11 o 3 o 

- ~ A X 2 Y + — X Li 2 (x) - 8 X ( 3 + - X 2 Li 2 (x) + — X vr 2 + — U X 2 - - X Y 3 

7 2 9 35 , x 15 2 , x / x 121 11 , 865 2 

-4 * 2 X 2 + - Li 2 (y) y - T vr 2 Li 2 (x) + 10 Li 3 (y) X - — U X - - X 4 - — X 2 



14 



+ Hx 3 --U 2 (y)Y 2 + -X :i Y + — Y X + -Y 2 X -- X 2 Y 2 + -tx 2 XY 
72 2 2yy! 6 72 3 8 4 



t 2 -s 2 



St 

2 



9 36 27 40 9 6 108 



22 



+2 U 4 (z) + 2 Li 4 (y) - y Li 3 (x) + 

41 ~ 



30377 ftT . , 440 rr ^ 

+ 2 Li 3 (x F + — c/y + 



155 



9 



1? 



Vy 



22 



-6 Cs ^ + - n 2 Y 2 - 2 X Li 3 (a;) -yIT + yl Li 2 (x) + 4 X Cs + X 2 Li 2 (x) 

+ ^X* 2 + 2 4uX 2 + 2XY 3 - 1 -n 2 X 2 + n 2 U 2 (x) + ^UY 2 -2 U 3 (y) X 
36 3 3 3 

93993 9 3 2 

+ ^XY--UXY + 2U 2 (y)XY-Y 4 --Y 3 -—Y 2 + -Y^ ^ " 
3 vy; 9 9 27 i 



if 



+ 1 



/2995 
432 



38 



509 



83 



83 



X + — vr 4 + — vr 2 + — Cs + 20 Li 4 (z) - — Li 3 (y) - 16 Li 4 (x) - 20 Li 4 (y) 



45 



48 



12 



Li 3 (x) + 10 Li 3 (x) y + ^ vr 2 y + 20 Li 3 (y) Y - 10 Cs Y + 2 tt 2 y 2 - 2 X Li 3 (x) 
12 6 

+ yX 2 y + gxLi 2 (x)+4XC3 + ^X 7 r 2 + 2Xy 3 + l 7 r 2 X 2 + ^Li 2 (y)y 



+4tt 2 U 2 {x) - 10Li 3 (y)X 



18 3 48 24 



^rx + -y 2 x + 2X 2 y 2 

24 12 



t 2 -s 2 



(3.34) 



11 



14 



119 



Un 2 + — X + 
6 9 180 



vr 4 + 24 tt 2 + 



173 



Cs- 



11 65 65 



299 

'IT 



55 



29 



1 



Li 3 (x) - 27Li 3 (:r) Y + — U Y - — vr 2 Y - 11 Cs Y + - A tx 2 Y 2 



+18XLi 3 (x) 



197 



X 2 Y 



3 
299 



6 

X Li 2 (x) + X Cs - 7 X 2 Li 2 (x) - y X vr 2 



6 6 

_ 11 [/ x 2 - — X Y 3 - - vr 2 X 2 - - vr 2 Li 2 (x) - — U Y 2 + 9 Li 3 (y) X - — U X 
6 6 4 4 v ; 3 v ; 6 

_ + !2 X 3 + 37 x3y _ m yA . 73 y2jf _ 55 

24 36 18 6 9 3 8 



+ 



/ 164771 121 



34 



- — U 2 + lOC/Ca /' ;; 2 

2592 9 S 9 



5941 



4 



1123 



■X- -7T 4 + 

432 3 144 



TT 



1529 
l36~ 



15 



- Yi? U + Y LU{z) + Y Li ^ } + 1 Li3 (a;) + if U Y ~ If Y ~ 38 Ca F 

A^Y 2 -23XU 3 (x) + ^X 2 Y-^XU 2 (x)+A2XC 3 + ^X 2 U 2 (x) 

_£1 x vr 2 + — [/ X 2 - — X F 3 + - vr 2 X 2 + — 7T 2 Li 2 (x) - — U Y 2 - 46 Li 3 (y) X 
72 3 6 12 4 v ; 3 vy; 

_l!l aX + « X 4 + 347 179 19 1169 59 

9 24 48 72 3 72 12 

+ E X 2 Y 2 + -7T 2 XY + -U XY + -Li 2 (y)XY + -Y A --Y 3 
8 4 3 2 vy; 12 18 

1169,. 5941 ^ 

+ — — Y 2 + -r-r-Y 



t 2 + s 2 ' 



72 216 
/ 11 TT o 207 v 

^7 



319 
~90~ 
121 

"lY 



1291 



121 



87 



^ 4 + T7T^ + ^rC3-^Li 4 (,) 



121 



144 



12 



Li 3 (2/) 



163 

+74 Li 4 (x) + ^ Li 4 (y) - Li 3 (x) - 26 Li 3 (z) F + — tt 2 F - 52 Li 3 (y) Y 



19 163 121 

+26 Cs F - — vr 2 Y 2 - 21 X Li 3 (rr) + — X 2 Y + — X Li 2 (x) - 34 X ( 3 



+-X 2 Li 2 (x) + -X7r 2 --f/X 2 + ixr 3 --7r 2 X 2 + — 
4 v ; 24 3 6 4 6 



Li 2 (y)F 



95 _2 T; , 0fiT: „ TTV 9^ 4 ,1291^2 433 w 17 



— 7r 2 Li 2 (a;) + 26Li 3 (y)X-imX--X 4 + m 



72 ^ + T Li 2 (y)^ 



+ 7x 3 r-^rx--r 2 x--x 2 F 2 + -7r 2 xr + -f/XF 

72 2 8 4 3 



17 > 
— Li 2 (y)XF 



t 2 -s 2 
st 



\ 12 3 360 



329 2 443 . 1502 TT 19139 

7T H 7T H C3 H - 

36 18 S 27 



324 



40 



+18 hu(z) + 18 Li 4 (y) + 18 Li 3 (x) F + vr 2 F - 18 Cs Y + 4 tt 2 F 2 - 18 X Li 3 (x) 



~ X 2 F + 18 X C 3 + 9 X 2 Li 2 (ar) + ^ X vr 2 + ^ tt 2 X 2 + 9 tt 2 U 2 (x) - 18 Li 3 (y) X 



11 

y 



+-X 4 -— X 2 + — X 3 + X 3 F + 7 F X + — X 2 F 2 + vr 2 X F 
2 36 3 2 



+ 18 Li 2 (y) X F - 7 F 2 - ^ F 



t 2 + s 2 



11 TT 2 121 v 119 4 157 2 71 . 43 T . . , 71 T . . , 

+ \ —Un 2 X tt 4 vr 2 C3 LiJz) H Li 3 (y) 

V 6 36 180 18 6 S 2 w 3 KyJ 

43 71 39 

+22 Li 4 (x) + — Li 4 (y) + — Li 3 (x) - 9 Li 3 (x) F - — vr 2 F - 18 Li 3 (y) F + 9 ( 3 Y 

17 3Q 71 Q 

-IL ^ Y 2 + 2 X Li 3 (x) - y X 2 F - ^ X Li 2 (x) - X C 3 - I X 2 Li 2 (x) 

* - 2 + H £/ X 2 - I X F 3 - H ^ X 2 - £ Li 2 („) F - \, 2 Li 2 (x) 



16 



+9Li 3 ( y )X + ^ X -gx 4 -^X 2 



^l X 3 --U 2 (y)Y 2 + -X 3 Y 
18 2 2Kyj 6 



lE Y X + -Y 2 X--X 2 Y 2 --7r 2 XY--UXY 
9 3 8 12 3 



+-Li 2 ( y )xr 



t 2 



-y 2 



t 2 M 2 



+ -Y 2 + 2YX -7T 2 -X' 



S 2 U 2 



(3.35) 



/ 22 TT 2 3133 v 203 4 401 2 1 . 83 T . , , 23 83 T . , , 

(-y f/ff+ i4T x -i8o^ + i8- ff -4 C3 + y Ll ^ ) + y + Y Ll ^ 

13 44 41 13 

+— Li 3 (x) + 41 Li 3 (x) Y + — 17 Y + - tt 2 Y + 9( 3 Y - -n 2 Y 2 -26X Li 3 (x) 

4 o y 12 

+yX 2 y-^XLi 2 ( a ;)+XC3 + |x 2 Li 2 ( a ;)-^Xvr 2 -| f /X 2 + ^Xy 3 

+| vr 2 X 2 + | ^ Li 2 (x) - f £/ y 2 - 1 1 Li 3 (y) X - f U X + | X 4 + ^ X 2 
209 X 3 _47 x3r _2^ rx _m F2x+ 125 x2y2 _33 ff2xF 



72 6 72 6 8 4 

+— t/ A" Y + — Li 2 ( s) XV- - F 4 + Hi y . + 259 y2 _ » \ 
3 2 2Vy; 4 9 72 72 J 

(rrr, 145 rr 2 743 v 113 4 1631 2 4 ^ 5297 rr 19 . . 

+ ( 5f/ C3 -^ f/ff+ l8- X - 240^ + W^- 9^ + ^ U + Y U ^ 
1Q 41 SQS 91 ^ 

+ Y LU(y) " ^592" + T Lis(;r) + 19 Lis {X) Y + 1 1 U Y + 2 ^ Y ~ G ^ 
+ A tt 2 Y 2 - 7 X Lia(x) + y X 2 y - ^ X Li 2 (x) - 2 X C 3 + ^ X 2 Li 2 (x) 

- ^ X vr 2 - H rj X 2 + I X y 3 + - A tt 2 X 2 + - A rr 2 Li 2 (x) - 1 1 U Y 2 + 5 Li 3 (y) X 
8 2 6 4 4 

11 7 a 97 9 17 o 4 „ 145 9 57 9 9 

— [/ x + — x 4 + — x 2 — x 3 --x 3 y y x - is y 2 x + — x 2 y 2 

2 12 48 8 3 24 8 

JI^xy + UUXY + -Li 2 ( y )xy - Ay 4 + 12 y 3 + —y 2 

12 2 vy; 12 24 



743 y 



24 



t 2 + s 2 

St 



11 „ 9 67 „ 14 4 1141 2 9 t 21 T . , s 9 



-y Li 4 (y) - Li 3 (x) + 5 Li 3 (x) Y + y vr 2 Y + 10 Li 3 (y) Y - 5 ( 3 Y + K 2 Y 2 

17 9 7 19 11 

-3 X Li 3 (x) + — X 2 y + - X Li 2 (x)+6I( 3 + -I 2 Li 2 (x) + — X tt 2 - — [/ X 2 
3 4 4 72 6 

+\XY 3 + 3 - t, 2 X 2 + | Li 2 ( B ) y + A i 2 Li 2 (x) - 5 Li 3 (») X+^UX+j X 1 

+^-X 2 - —X s + - U,(y) Y- 2 -lx>Y-™YX-* Y 2 X + 5 A" 2 F 2 
144 72 2 vy; 3 72 12 8 



17 



1 11 1 \ It 2 - s 2 " 

--^XY+-UXY-- U 2 (y) X Y'\ \—- 



( 21 47 4 47 2 11 11 , . 11 t ' / \ 11 T , x 

+ ~ 16 ~ ^80 16 T T W 2 ) + y Li 3 (y) + 6 Li 4 (z) + — Li 4 (y) 

+ H Li 3 (x) - 3 Li 3 (x) y - | vr 2 y - 6 Li 3 (y) Y + 3 Cs Y - I vr 2 y 2 - | X 2 y 

" y X Li 2 (x) - 3 X Cs - ^ X 2 Li 2 (x) + l l X vr 2 - I X y 3 - ^ vr 2 X 2 - H Li 2 (y) y 

- 2 Li 2 (x) + 3 Li 3 (y) X - 1 X 4 - ^ X 2 - | X 3 + i Li 2 (y) y 2 + 1 x 3 y 

+ t Y x ~ \ y2 x - \ x2y2 + r 2 x Y ~ \ Li2 ^ x Y ) [~u^ 

+ y 2 ^ + (y 2 - 2 yx + 7 r 2 + x 2 )-^ (3.36) 



1 29 v 3 4 25 
- H X + -7r 4 H 

2 4 5 6 



7T 2 + 15 C 3 - 10Li 4 (V) - 10Li 4 (y) - 19Li 3 (x) 



+12 Li 3 (x) y + - vr 2 y + 12 C 3 Y + - vr 2 y 2 + 4 X Li 3 (x) + 13 X 2 y + 19 X Li 2 (x) 
3 3 

-16 X Cs + X 2 Li 2 (x) - - X vr 2 + 2 X y 3 - - vr 2 X 2 + tt 2 Li 2 (ar) + 20 Li 3 (y) X 

6 3 

-x 4 + 4 x 2 - — x 3 + 4 x 3 y - y x - 9 y 2 x + - x 2 y 2 - - tt 2 x y 

6 2 3 



+2 Li 2 (y) x y - y 4 + 6 y 3 + y 2 

+ f -6Z7C3 



29 



y 



1 rr 2 4 7 v 7 4 3 2 23 A 3 TT 11 T . . . 
v 2^ + 16 X+ lT "IT + T C3 -8 f/ -Y Ll ^ ) 

-y LU(y) - f Li3(x) - ^ + 7Li 3 (x) y + ^ 2 y + 7C,3 Y + A vr 2 y 2 

+2 X Lia(x) + ^ X 2 y + H x Li 2 (x) - 9 X ( 3 + \ X 2 Li 2 (x) + ^ x tt 2 + ^ X y 3 

-A. 2 X 2 + ^ 2 Li 2 (a;) + llLi3(y)X-^X 4 -^X 2 -^X 3 + ^X 3 y + |yX 

-5y 2 x + A x2r2 _^ 2xr + 3 Li2(y)xr _^ r4+ w r 3_^ y2 



t 2 + s 2 

St 



18 



/ 31 „ 103 , 23 , 19 , 11 i 19. . . . 

+ ( - T6 X + 180 * - 16 * + T Ca + T Lu(z) - T Ul(y) - 14 Lulx) 

1119 1 
-— Li 4 (y) - — U 3 (x) + 3 Li 3 (x) Y + 2 ir 2 Y + 6 Li 3 (y) Y - 3 Cs ^ + - tt 2 F 2 

19 17 1 

+4 X Li 3 (x) + 2 X 2 y + — X Li 2 (x) + 7X( 3 + ^ 2 Li 2 (x) + — X tt 2 + - X Y 3 

+ 3 - n 2 X 2 + ^ Li 2 (y) Y-\, 2 Li 2 (a;) - 3 Li 3 (y) X - 1 X 4 - § X 2 + ^X 3 

-iLi 2 ( y )y 2 + ix 3 y + |yx + H r2x+ 5 x2y 2_i9 7r2xr 



+ l -U 2 {y)XY 



t 2 -s 2 

si 



+ (2X 2 -yx + -7r 2 + y 2 



t 2 + s 2 



+ (^X 2 -3Y X + ^tt 2 -3X 



t 2 -s 2 
■u 2 



y 2 



t% 2 



+ 



y 2 + 2yx-7r 2 -x 2 



S 2 U 2 



(3.37) 



X tt 4 tt 2 -4Cs + -C/-4Li 3 (x) -2Li 3 (a;)y + — 17 Y 



/91 11 4 25 , . . 2„ ... , , , 10 

X 7T 

V 18 180 18 

2 1 37 2 
+-7T 2 y-2C 3 y + -7r 2 y 2 -- + 2X 2 y + 4XLi 2 (:r) + 2XC 3 + -X7r 2 

o O lo o 

-ixy 3 -2Li 3 (y)x--[/x-2X 2 + 4yx-2y 2 x-ix 2 y 2 + -7r 2 xy 

3 3 2 3 

+ ±Y 4 + -Y 3 -4Y 2 -—y) 
6 3 9 



44 „„ 11 „ o 5„ 11 „ 133 2 13 , 185 rr 2 T ... 

7T ( 3 H [/ Li 3 (x) 

6 54 3 v ; 



+ 1 -— U 2 + — Un 2 + -X + — tt 4 + 
V 9 36 9 90 108 



+4Li 3 (a;)y + ^f/y + ^7r 2 y + 4C3y-^7r 2 y 2 -ix 2 y + ^XLi 2 (a;) 
9 36 3 4 3 



1 



22 



4 X Cs + — X tt 2 + - [/ X 2 + - X Y 3 + - U Y 2 + 4 Li 3 (y) X [/ X 



18 



9 



- ' v3 i V V i v2 "V i V"2 v2 _ 7r 2 J^y_ — [/ X Y 



4 
3 



2 
3 



y 4 



6 



y J 



— y^ + 

18 216 



1307 10 



9 



y 



t 2 + g 2 

si 



+ (- ^vr 2 -^X + ^7r 4 + ^7r 2 -^C3 + 8Li 4 (^) + i Li 3 (y) - 8Li 4 (x) 



-8 Li 4 (y) + | Li 3 (x) + 4 Li 3 (x) y + - A tt 2 y + 8 Li 3 (y) Y-4( 3 Y + n 2 Y 2 
+ ^X 2 y-^XLi 2 (x)+4XCs 



18 



Xtt 2 



X o 2 o 2 o f 

-f/x 2 + -xy 3 + -7r 2 x^ 



4 4 22 1 79 7 

-- Li 2 (y) y + - tt 2 Li 2 (x) - 4 Lia(y) X + y f/X + -X 4 + -X 2 --X 3 



19 



-X 6 Y - — YX - — Y 2 X + X 2 Y 2 + -UXY ) 
3 18 12 3 



+ 1 
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4. One-loop self-interference 

In this section, we present the self- interference of the one- loop amplitude of Eq. (2.17) 
in terms of the one-loop box in D = 6 — 2e dimensions and the one-loop bubble in 
D = A — 2e dimensions. The e expansion of these integrals is given in the Appendix 
and can be inserted directly into the following expressions. We can write 

C 8(lxl) (s, t, u) = -V s {s, t, u) T 2 (s, t, u) sHu + 2 Re {V c {s, t,u)} + T T {s, t, u), (4.1) 

where the self-interference of the singular terms of the amplitude is 
VN, 



V s (s,t,u) 



■\V c (s,t,u) 



+ j2 Re {V c (s, t, U y [P A (s, t, u) + TZ A (s, t, u) + (t <-> u)} } 



V 2 



AN 
V 

"an 



{\v a (s, t, u) + n A (s, t, u) | 2 + {t <-> U)} 

2 Re { [v A (s, t, u)t + K A (s, t, w) f ] [V A (s, u, t) + K A (s, u, t)}} 



(4.2) 



and the interference of the singular terms with the one-loop amplitude is 
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u 
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+ 4 
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V 2 



[(V A (s, t, w) f + TZ A (s, t, u)^ uC 2 {s, t, u) + (t<-+ u) 
[(V A (s, u, + K A (s, u, t) f ) uC 2 (s, t, + u)] ] , (4.3) 



AN 
V_ 

~AN 

where the functions C 2 and £3 are defined in Eqs. (3.13) and (3.14) respectively. In 
Eqs. (4.2) and (4.3) 
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originates from ultraviolet terms which are rendered finite after renormalisation in 
the MS scheme. 

Finally, the finite contribution is given by 

VN V V 2 V 
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with 



T e (t, s) = Bub (t) - (1 + e)Bub (s) . 



(4.12) 



Only the first term in the expansions of the box integrals are required in Eqs. (4.8)- 
(4.11), where we have systematically discarded contributions of 0(e). The bubble 
integrals must be expanded through to O (e°). 

In this section, we have isolated the infrared divergences of the one-loop ampli- 
tude into the terms Va and Vc- Both functions depend on one-loop bubble integrals 
and diverge as 1/e 2 . The separation of the singular terms becomes evident when we 
complete our basis of one-loop master integrals with the finite box in D = 6 — 2e. 
In the context of one-loop integrals, this form for the divergences arises naturally. 
The singular behaviour of the one-loop amplitude can also be predicted applying the 
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formalism of [25] which yields 



V AiC (s,t,u) = --A(e,s,u,t) (4.13) 

and 

V c ,c(s, t, u) = B(e, s, t, u) - -B(e, s, u, t), (4.14) 

su st 

where A(e,s,t,u) and B(e,s,t,u) are defined in Eqs. (3.3) and (3.4). The singular 
structure of Va,c an d ^Pc,c around e = precisely matches that given in Eq. (4.4) 
and (4.5) and the two formalisms differ only in the finite remainder. In order to 
make the agreement explicit one could replace Va{Vc) with Va,c{^Pc,c) m Eqs. (4.2) 
and (4.3) with appropriate modifications due to the finite differences V A — Va,c an d 
V c ~ V C)C in E qs- ( 4 -3) and (4.7). 



5. Summary 

In this paper we presented the O (aj) QCD corrections to the 2—^2 scattering 
processes qq — > gg, gg — > qq and the associated crossed processes qg — > qg and qg — > 
qg in the high energy limit, where the quark masses can be ignored. We computed 
renormalised analytic expressions for the interference of the tree-level diagrams with 
the two-loop ones and for the self-interference of one-loop graphs in the MS scheme. 
Throughout we employed conventional dimensional regularisation. 

The renormalised matrix elements are infrared divergent and contain poles down 
to O (1/e 4 ). The singularity structure of one- and two-loop diagrams has been thor- 
oughly studied by Catani [25] who provided a procedure for predicting the infrared 
behaviour of renormalised amplitudes. The anticipated pole structure agrees exactly 
with that obtained by direct Feynman diagram evaluation. In fact Catani's method 
does not determine the 1/e poles exactly, but expects that the remaining unpredicted 
1/e poles are non- logarithmic and proportional to constants (colour factors, it 2 and 
( 3 ). We find that this is indeed the case, and the constant is given in Eq. (3.19). 
This provides a very strong check on the reliability of our results. Similarly, the 
infrared divergent structure of the squared one-loop diagrams we found by direct 
evaluation agrees with the expected pole structure. 

The pole structure of the two-loop contribution is given in Eq. (3.2) while ex- 
pressions for the finite parts in the s-channel and ^-channels according to the colour 
decomposition of Eq. (3.24) are given in Sees. 3.2.1 and 3.2.2 respectively. Similarly, 
the infrared divergent one-loop contributions along with the remaining finite parts 
are detailed in Sec. 4. The expressions for the two-loop pole structure and the sin- 
gular and finite parts of the self-interference of one-loop graphs are analytic and are 
given in terms of the one-loop bubble graph and one-loop box graph in D = 6 — 2e 
dimensions. To evaluate these formulae in the appropriate physical region requires 
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the insertion of the series expansion of the one- loop graphs around e = 0. These 
expansions are given in Appendix A. 

In this paper, we have concentrated only on QCD processes. However, the QED 
processes e + e~ — > 77 and Compton scattering e~7 — > e~ r y as well as 77 — > e~e + are 
also of interest. We note that the expressions given here are related to those for the 
corresponding QED processes in the limit where the electron mass can be ignored 
by the alteration of various colour factors. We expect to address this in a separate 
article. 

The results presented here, together with those previously computed for quark- 
quark scattering [22, 23, 24] are necessary ingredients for the next-to-next-to-leading 
order predictions for jet cross sections in hadron-hadron collisions. On their own, 
they are insufficient to make physical predictions and much work remains to be 
done. First, at the matrix element level, similar expressions to those presented here 
for gluon-gluon scattering are needed. Given the recent progress in the field, we 
anticipate that this problem will soon be solved. Second, a systematic procedure 
for analytically canceling the infrared divergences between the tree-level 2^4, the 
one-loop 2-^3 and the 2-^2 processes needs to be established for semi-inclusive jet 
cross sections. Again, recent progress in determining the singular limits of tree-level 
matrix elements when two particles are unresolved [36, 37] and the soft and collinear 
limits of one-loop amplitudes [38, 39], together with the analytic cancellation of the 
infrared singularities in the somewhat simpler case of e + e~ — > photon + jet at next- 
to- leading order [40], suggest that the technical problems will soon be solved for 
generic 2 — > 2 scattering processes. There are additional problems due to initial 
state radiation. However, the recent steps taken towards the determination of the 
three-loop splitting functions [3, 4, 5] are also promising. 

Third, a numerical implementation of the various contributions must be devel- 
oped. The next-to-leading order programs for three jet production that have already 
been written provide a first step in this direction [41, 42]. We are confident that 
the problem of the numerical cancellation of residual infrared divergences will soon 
be addressed thereby enabling the construction of numerical programs to provide 
next-to-next-to-leading order QCD estimates of jet production in hadron collisions. 
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A. One-loop master integrals 



In this appendix, we list the expansions for the one-loop box integrals in D = 6 — 2e. 
We remain in the physical region s > 0, u,t < 0, and write coefficients in terms of 
logarithms and polylogarithms that are real in this domain. More precisely, we use 
the notation of Eqs. (3.22) and (3.23) to define the arguments of the logarithms and 
polylogarithms. The polylogarithms are defined as in Eq. (3.21). 
We find that the box integrals have the expansion, 



Box 6 (w,t) = 



e^r(l + e)r(l -e) 2 //i 2 ' 



2sY (1 - 2e) (1 - 2e) 
+2e 



U 3 (x)-XU 2 (x)-^X 3 -^X 



(X-YY + tt' 

2 



-2e 2 



Li 4 (a;) + YLi 3 (x) - -X 2 Li 2 (x) - -X 4 - -X S Y + -X 2 Y 2 



-^X 2 - n -XY-^ 
4 3 45 



+ (u <-> t) 1 + O (e 3 ) , (A.I) 



and 

Box 6 (s,t) 



e"T(l + e)r(l-e) 2 
2ur(l-2e)(l-2e) \ u 



X 2 + 2inX) 



+e 



(-2Li 3 (x) + 2XLi 2 (x) - jjX 3 + 2YX 2 - tt 2 X + 2( 3 ) 



+m 2Li 2 (:r) +4FX -X 2 - 



7T 



2Li 4 (z) + 2Li 4 (y) - 2YU 3 (x) - 2XLi 3 (y) + (2XY - X 2 - 7r 2 )Li 2 (:r) 



+ \x A - |x 3 F + -X 2 F 2 + Jtt 2 X 2 - 2tt 2 XF + 2FCs + ^ 



2 
3' 



1 

6 : 



+in( K -2U z (x) - 2Li 3 (y) + 2FLi 2 (:r) + ^X 3 - 2X 2 F + 3XF 2 

■ O (e 3 ) . (A.2) 



7T 



Box 6 (s,w) is obtained from Eq. (A.2) by exchanging u and t. 

Finally, the one-loop bubble integral in D = 4 — 2e dimensions is given by 



Bub(s) = 



e^r(l + e)r(l-6) z / /i 2 
r (2 - 2e) e V s 



(A.3) 
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